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Abstract 


The  purpose  of  this  research  was  to  search  for  evidence  of  low  temperature 
annealing  from  neutron  irradiated  4H-silicon  carbide.  No  features  suggesting  annealing 
were  found  below  a  temperature  of  340K.  Temperature  dependant  Hall  effect 
measurements  were  taken  over  a  range  of  lOOK  to  340K  recording  resistivity,  carrier 
densities,  and  mobility.  Resistivity  was  noted  to  increase  with  irradiation,  and  carrier 
densities  appeared  to  decrease,  while  mobility  appeared  minimally  affected  by  neutron 
irradiation.  This  suggests  the  creation  of  active  acceptor  defects  decreasing  carrier 
concentrations.  N-type  samples  measured  were  5mm  x  5mm  square  with  Nickel  contacts, 
and  irradiated  to  lO'^  and  lO'^  n-cm/cm^  of  IMeV  equivalent  neutron  fluence. 
Suggestions  for  continuing  research  include  using  a  probe  station  instead  of  wire 
connections  to  samples,  use  a  large  source  current  to  minimize  variance,  and  minimize 
cadmium  shielding  to  reduce  negative  reactivity. 
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HALL  MEASUREMENTS  OF  SURFACE  DEFECTS  IN 
NEUTRON  IRRADIATED  SILICON  CARBIDE 

I.  Introduction 

The  Promise  of  Silicon  Carbide 

Just  as  silicon  revolutionized  electronics  in  the  twentieth  century,  silicon  carbide 
will  change  how  we  think  about  electronics  in  the  twenty-first  century.  Silicon  carbide 
(SiC)  is  a  wide  band-gap  semiconductor,  the  material  of  choice  for  extreme  environment 
electronics.  Silicon  carbide  is  attractive  for  use  in  high-temperature,  high-power,  high 
frequency,  and  highly  radioactive  environments.  The  United  States  Air  Force  is 
researching  methods  of  exploiting  the  material  to  improve  access,  reliability,  and  lethality 
of  military  assets.  This  thesis  proposes  to  look  at  the  effects  of  neutron  damage  to  silicon 
carbide. 

Silicon  carbide  is  attractive  for  use  in  high-temperature,  high-power,  high- 
frequency,  and  highly  radioactive  environments.  Traditional  silicon  has  a  lifetime  of  10 
years  in  an  environment  of  100°  C,  and  is  unable  to  operate  above  250°  C  [36,4].  SiC  is 
expected  to  survive  for  a  similar  lifetime  in  an  environment  of  500°  C  and  higher  [36]. 
Additionally,  silicon  carbide’s  large  breakdown  electric  field  suggests  excellent 
suitability  for  high-frequency  applications  [48].  McLean  note  survivability  of  SiC  devices 

up  to  a  neutron  fluence  of  lO'^ - ^  [36].  With  respect  to  primary  damage  caused  by 

cm 

neutron  irradiation,  SiC  is  superior  to  silicon  [26]. 

To  the  military,  silicon  carbide  means  easier  maintenance,  faster  aircraft,  and 
improved  survivability  for  space-based  platforms.  Silicon  carbide’s  high  thermal 
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conductivity  allows  for  device  plaeement  in  engines  decreasing  troubleshooting  time.  SiC 
deviees  will  reduee  eooling  requirements,  making  platforms  lighter,  increasing  speed  and 
deereasing  mass  budgets.  SiC  promises  to  “provide  high  power  solid-state  radars  for 
future  taetieal  systems”  [9]. 

The  USAF  understands  the  potential  of  SiC  devices  and  has  devoted  appropriate 
resources  to  developing  the  material.  From  1992  to  1996  the  USAF  established  the 
eonsortium  for  the  development  of  silieon  earbide  for  eleetronies  applieations,  a  group  to 
study  SiC  supporting  researeh  at  US  Universities  [9].  The  Air  Foree  Researeh  Laboratory 
(AFRL)  eurrently  studies  advaneed  eontaets  on  SiC  as  the  promise  of  SiC  based  deviees 
is  hindered  hy  the  limitations  of  deviee  eontaets.  Multiple  theses  eompleted  at  the  Air 
Foree  Institute  of  Teehnology  eoneentrate  on  eharacterizing  the  properties  of  SiC. 

The  interest  in  silieon  carbide  surfaces  in  partieular  is  the  result  of  device 
vulnerabilities  at  junetions.  For  junetion  deviees,  ionizing  radiation  ineident  on  active 
junctions  results  in  an  increase  of  earriers.  For  metal  oxide  semiconductor  (MOS) 
devices,  ionizing  radiation  results  in  eharge  trapping  at  the  oxide-semieonduetor  junetion 
[17:134].  Non-ionizing  radiation,  caused  by  interaetions  with  heavy  particles  leads  to 
displaeement  damage  such  as  point  defects  in  the  erystal  structure.  These  defects  in  the 
lattice  eause  eharge  earriers  to  seatter,  deereasing  mobility  and  increasing  resistivity. 

Objective 

This  experiment  attempts  to  examine  the  aimealing  behavior  of  the  surface  of 
neutron  damaged  4H-SiC  over  a  temperature  range  from  100  K  to  350  K. 

Charaeterization  of  the  post-irradiated  surface’s  eleetrieal  properties  will  aid  in 
identifying  modes  of  failure.  By  eomparing  temperature  dependency  of  mobility  in 
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irradiated  and  non-irradiated  samples,  displaeement  damage  ean  be  identified  by 
observing  reeovery  of  resistivity,  mobility,  or  earrier  eoneentration. 

Paper  Organization 

The  following  thesis  is  divided  into  supporting  theory,  experimental  proeedure, 
results,  and  eonelusion.  Chapter  II  introduces  the  theory  behind  models  of  conduction, 
the  crystal  structures  of  SiC,  and  a  primer  on  Hall  Effect  measurements.  Chapter  III 
discusses  the  experiment  including  equipment,  the  planned  experiment,  and  the  changes 
that  were  required.  The  results  from  the  experiment  are  presented  and  analyzed  in 
Chapter  IV.  Chapter  V  extracts  the  most  fundamental  results,  and  proposes  courses  of 
action  based  upon  knowledge  gained. 
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II.  Theory 

Devices  made  from  silicon  carbide  are  more  radiation  tolerant  than  traditional 
silicon  devices,  making  SiC  attractive  for  use  in  the  space  and  nuclear  reactor 
environments.  Both  environments  are  filled  with  higher  energy  particles  than  normally 
found  on  the  earth’s  surface.  The  next  section  describes  how  heavy,  neutral  particles 
induce  damage  to  SiC. 


Heavy  Particle  Interactions  with  Matter 

Heavy  Particles,  such  as  neutrons,  protons,  and  alpha  particles,  can  impart  enough 
energy  through  collisions  to  remove  an  atom  from  its  position  in  a  crystal  lattice,  creating 
a  point  defect.  These  defects  are  often  called  Frenkel  pairs,  which  are  a  pair  of  simple 
defects  -  a  vacancy  and  an  interstitial  [17:73].  Frenkel  defects  can  create  additional 
energy  levels  within  the  band  gap,  resulting  in  trap  sites  [46: 11].  Heavy  particles  transfer 
only  a  fraction  of  their  kinetic  energy  in  a  single  collision  and  below  certain  energies,  a 
particle  carmot  impart  enough  energy  to  cause  a  displacement.  Equation  1  describes  the 
relationship  between  the  mass  of  both  particles,  the  energy  of  the  incident  particle,  and 
the  maximum  fraction  of  energy  imparted  in  the  collision.  The  first  atom  hit  in  the  lattice 
by  a  particle  is  often  referred  to  as  the  primary  knock-on  atom  (PKA).  The  PKA  often 
absorbs  enough  energy  to  create  many  secondary  displacements  in  the  lattice.  Equation  1 
also  describes  the  maximum  fraction  of  energy  imparted  to  a  secondary  atom  by  the 
PKA. 


=  - - —  [54:213] 

En  {m  +  My 


(1) 
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where,  m,  M  Masses  of  the  particles  involved,  and 

Qmax/En  ~  Fraction  of  Energy  imparted  to  PKA. 

Given  a  neutron  with  kinetic  energy  E,  a  maximum  of  0.1 33E  is  imparted  to  a 
silicon  atom,  and  0.284E  to  a  carbon  atom  as  calculated  through  use  of  (1).  The  threshold 
energy  for  an  incident  neutron,  below  which  displacement  is  unlikely,  is  10  keV  [17:73]. 
While  a  neutron  creates  only  a  small  number  of  PKAs,  each  PKA  has  the  potential  to 
create  thousands  of  secondary  defects  within  the  lattice. 

The  Kinchin-Pease  model  of  defect  formation  states  that  the  number  of  defects  is 
proportional  to  the  energy  of  the  PKA  [17].  Also,  the  number  of  Frenkel  pairs  is 
proportional  to  the  energy  given  to  the  PKA  [17:73].  Using  (1)  again,  we  find  that  an 
atom  can  impart  up  to  all  of  its  energy  to  another  atom  of  the  same  kind.  Since,  an  atom 
can  impart  all  of  its  energy,  a  40  eV  PKA  can  produce  up  to  1  defect.  By  the  Kinchin- 
Pease  model,  that  means  a  133  keV  PKA  can  create  more  than  3200  defects  in  SiC. 

A  particle  only  imparts  the  maximum  amount  of  energy  during  a  head-on 
collision.  However,  it  is  more  likely  that  a  collision  will  occur  at  another  angle  of  attack. 
The  energy  transferred  to  an  atom  by  a  neutron  is  more  generally  described  by  equation 

2. 


where. 


(2) 


E  Energy  of  neutron, 

A  Atomic  Mass  of  Target  Atom,  and 

0  Incident  angle  of  attack. 
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Since  the  ineident  angle  of  attaek  for  a  partieular  collision  is  unknown,  the 
average  energy  imparted  is  used  over  a  large  quantity  of  eollisions.  The  average  energy 
transferred  ( )in  a  eollision  is  aceepted  to  be  half  of  the  maximum  energy  as  denoted 
in  equation  2.  It  is  this  energy  that  is  the  quantity  of  interest  for  dose  caleulations. 


Qavc 


(3) 


Electrons  and  gamma  rays  often  do  not  transfer  enough  energy  to  create  a  Frenkel 
pair  owing  to  their  lower  momentum,  and  because  of  this  a  neutron  displaees  a  thousand 
times  more  atoms  than  an  equivalent  energetie  eleetron  [17:77].  Therefore  the  effects  of 
electrons  and  gamma  rays  will  be  ignored  in  this  study. 

The  dose  to  silieon  carbide  for  the  neutron  flux  of  the  reaetor  ean  be  estimated  by 
ealculating  the  “first-eollision”  dose  to  the  SiC  for  all  neutrons  with  E  >  10  keV  [54:371, 
17:73].  The  first-eollision  dose  sets  a  lower  bound  for  the  dose  because  it  assumes  that 
the  neutron  escapes  the  system  before  interaeting  more  than  once.  This  is  a  good  estimate 
for  systems  where  the  thickness  is  much  less  than  the  mean  free  path  length  of  the 
neutron.  While  reactors  have  different  neutron  energy  distributions,  comparisons  to  other 
radiation  environments  ean  be  made  by  normalizing  the  fluence  to  a  1  MeV  equivalent 
fluenee.  An  example  of  finding  IMeV  equivalent  fluence  is  supplied  in  Appendix  B. 

Silicon  carbide  ean  survive  a  higher  neutron  dose  than  silicon  due  to  stronger 
moleeular  bondings.  The  amount  of  energy  required  to  free  an  atom  from  the  lattice 
differs  based  on  the  atom  and  the  geometry  of  its  bonds.  It  takes  about  40eV  of  energy  to 
eause  a  displacement  in  SiC.  Table  1  exhibits  the  amount  of  energy  required  to  free  a 
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silicon  or  carbon  atom  in  specific  planes  of  hexagonal  silicon  carbide.  For  comparison,  a 
silicon  atom  in  a  silicon  lattice  requires  only  13-20  eV  to  be  displaced  [26]. 

Table  1:  Displacement  Energy  for  Si  &  C  atoms  in  the  different  crystal  plans  of  SiC. 

Compiled  from  [45]. 


Plane 

Silicon  (eV) 

Carbon  (eV) 

001 

30  -  35 

35  -  40 

111 

35 

20  -  25 

110 

80  -  85 

30 

Defects  in  crystal  structures  can  be  removed  through  annealing.  The  two  basic 
forms  of  annealing  are  isochronal  and  isothermal  [17:80-84].  Isochronal  annealing 
studies  a  range  of  temperatures  to  find  the  specific  temperatures  at  which  defects  become 
mobile,  returning  to  a  normal  lattice  position.  Isothermal  annealing  studies  the  time 
dependent  change  in  properties  at  a  specific  temperature.  For  this  study,  in  order  to 
identify  if  annealing  is  occurring  below  250°  C,  an  isochronal  annealing  study  was 
accomplished.  In  traditional  silicon,  Holmes-Siedel  notes  that  the  only  annealing  that 
occurs  at  temperatures  below  200°  C  is  at  the  semiconductor’s  surface  [17:83].  This  data 
should  lend  to  support  the  possibility  of  similar  annealing  in  silicon  carbide. 

Kronig-Penney  Model  of  Conduction 

Columbic  force  is  the  driver  that  shapes  crystal  structures.  The  Kronig-Penny 
model  crudely  models  the  Coulombic  potential  in  a  lattice  as  a  periodic  square  wave 
[35:321].  This  one-electron  approximation  assumes  that  the  random  interactions  of  other 
electrons  in  the  lattice  sum  to  zero,  so  that  the  periodic  potential  is  all  that  remains. 

Others  have  simplified  this  model  to  represent  the  potential  by  periodic  Dirac  delta 
functions  [14:200].  In  either  case,  it  is  the  periodicity  of  the  potential  yielding  useful 
information.  The  solution  of  the  Schrodinger  equation  with  these  periodic  potentials  leads 
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to  bands  of  allowed  and  forbidden  energy  states.  While  not  a  precise  approximation  of 
the  potentials  felt  in  a  crystal,  the  Kronig-Penney  model  has  shown  itself  to  be  a  useful 
tool.  The  Kronig-Penney  models  states  that  if  the  periodic  structure  is  interrupted  with  a 
defect,  the  electrical  properties  will  change. 

A  charge  carrier  traveling  through  a  crystal  with  no  defects  would  continuously 
accelerate  under  the  influence  of  an  electric  field.  However,  this  behavior  is  not  observed. 
Instead  crystals  have  many  types  of  defects  that  act  as  scattering  points  for  carriers. 
Defects  create  allowable  states  within  the  bandgap,  trapping  minority  carriers  from 
conducting.  As  defects  in  a  crystal  increase,  conductivity  decreases.  The  diversity  of  SiC 
polytypes  increases  the  difficulty  of  describing  the  material  in  general  as  each  crystal 
structure  leads  to  slightly  different  electrical  properties. 

Silicon  Carbide 

Silicon  carbide  has  been  the  subject  of  study  for  many  decades.  The  first  studies 
were  done  at  the  beginning  of  the  twentieth  century,  but  getting  high  purity  crystals  was  a 
problem  for  a  long  time.  In  the  1980’s  research  quality  6H-SiC  became  available 
commercially,  followed  by  4H-SiC  in  the  1990’s.  Due  to  availability,  most  research  has 
concentrated  on  these  two  polytypes,  preferring  4H  owing  to  its  larger  bandgap. 

Crystal  Structures 

Silicon  carbide  has  many  different  polytypes.  Polytypism  is  a  one-dimensional 
form  of  polymorphism  where  the  closest  packed  layers  are  similar,  and  crystals  vary  in 
the  pattern  of  stacking  layers  [48].  While  nearly  200  polytypes  of  SiC  have  been 
identified,  only  a  few  have  the  attention  of  researchers. 
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Silicon  carbide  lattice  structures  fall  into  three  main  types;  Cubie  (C),  Hexagonal 
(H),  and  Rhombic  (R).  The  stacking  sequence  of  a  polytype  is  usually  described  in 
Ramsdell  or  ABC  notation,  where  each  layer  arbitrarily  identifies  a  similar  layer.  The 
number  preceding  the  letter  C,  H,  or  R  describes  the  number  of  stacks  in  the  pattern 
followed  by  a  letter  deseribing  the  erystal  structure  [42:138].  Only  one  cubic  structure 
exists  for  SiC,  and  it  is  referred  to  as  both  p-SiC  and  3C-SiC.  (3-SiC  possesses  the 
common  zinc-blende  lattiee  structure  with  the  stacking  pattern  of  ABC.  The  two  most 
common  hexagonal  types  are  6H  and  4H  SiC.  6H-SiC  has  an  ABCACB  stacking  pattern, 
and  4H-SiC  has  an  ABAC  stacking  pattern.  Hexagonal  and  rhombic  SiC  are  sometimes 
colleetively  referred  to  as  a-SiC. 


Figure  1:  Structure  of  Common  SiC  Polytypes.  Source  [42:148] 

Each  crystal  lattice  of  SiC  has  a  unique  bandgap  that  ranges  from  2  to  4eV.  4H- 
SiC  is  the  preferred  crystal  for  researchers  because  its  bandgap  is  larger  than  6H-SiC. 


9 


Activation  Concerns 

Silicon  carbide  interests  nuclear  researehers  in  partieular  due  to  its  low  neutron 
absorption  cross  section.  Silicon  and  carbon  both  possess  small  neutron  absorption  cross 
sections.  Additionally,  most  silicon  and  carbon  atoms  must  absorb  2  or  3  neutrons  before 
beeoming  an  active  isotope.  When  silicon  activates  as  Si-31,  it  decays  quickly  to  a  stable 
isotope  of  phosphorous.  This  charaeteristie  ean  be  used  to  dope  SiC  with  silicon  through 
a  process  known  as  neutron  transmutation  doping  (NTD).  The  activation  of  SiC  comes 
mostly  from  C-14,  whieh  has  a  long  half-life  of  5715  years.  For  the  reason  of  low 
aetivation,  SiC  composites  are  being  developing  as  struetural  materials  for  possible 
nuelear  fusion  reactors.  After  30  years  of  operation  in  a  fusion  reactor,  SiC  could  be 
safely  disearded  by  shallow  burial. 

For  this  researeh,  SiC’s  low  activation  potential  means  that  the  main  activation 
concern  stems  from  dopants,  ohmie  contacts,  and  wire  bonds.  Nitrogen  and  aluminum  are 
the  dopants  for  n-  and  p-type  SiC  respectively.  Nitrogen  possesses  a  small  absorption 
eross  section,  decaying  quickly  to  oxygen,  and  is  of  little  coneem.  Aluminum  also  decays 
quickly,  to  silicon.  Ohmic  materials  selected  were  niekel,  aluminum,  and  titanium  for 
reasons  eonnected  to  ability  to  form  a  good  contact. 

Of  the  mentioned  ohmie  contaet  materials,  none  pose  a  potential  of  eausing  an 
aetivation  problem.  While  ^  of  nickel  will  transmute  to  Ni-59  with  a  half-life  of  7600 

years,  it  deeays  by  eleetron  eapture  and  only  emits  a  neutrino.  While  A1  activates  more 
frequently,  it  deeays  quiekly  beeoming  insignifioant  within  an  hour. 

While  the  ereation  of  low  resistance  ohmic  contacts  drove  the  selection  of  contact 
material,  aetivation  was  the  driving  force  in  selecting  wire  material.  Potential  materials 
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were  compared  using  neutron  absorption  cross  sections,  and  assuming  transmutation  to 
the  n+1  isotope.  By  making  this  assumption,  the  (n,p)  reaction  was  ignored.  Ignoring 
other  nuclear  reactions  underestimates  total  activation.  Table  2  contains  common 
elements  used  in  electronics,  with  cadmium  listed  since  it  is  often  used  to  absorb  lower 
energy  neutrons.  Elements  with  “Mostly  Stable”  listed  under  half-life  need  to  capture  two 
or  more  neutrons  before  they  transmute  into  a  radioactive  isotope.  The  daughter  products 
of  all  the  listed  elements  are  stable.  Silicon  and  carbon  both  rarely  transmute  to  an  active 
radioisotope.  When  they  become  active,  only  C-14  maintains  a  high  activity  for  a 
significant  time,  as  Si-31  has  a  half-life  of  2.6  hours. 


Table  2:  Thermal  Cross  sections  and  Half-Lives  of  elements  common  in  electronics. 

Compiled  from  [2:68,  43] 


Element 

Thermal  Crossection  [barns] 

Half-Life  (after  n-absorbtion) 

Cadmium 

2870 

Mostly  Stable 

Silicon 

0.14 

Mostly  Stable 

Carbon 

0.0035 

Mostly  Stable 

Nickel 

4.1 

Mostly  Stable 

Nitrogen 

1.67 

Mostly  Stable 

Aluminum 

0.217 

2.25  min 

Titanium 

5.14 

Mostly  Stable 

Gold 

98.7 

2.7  days 

Indium 

87 

54  min 

Gold  and  indium  are  elements  common  for  wire  bonding.  The  combination  of  a 
relatively  large  thermal  absorption  cross  section  and  a  2.7  day  half-life  make  gold 
undesirable  for  use  in  this  project.  Since  a  wire  bonder  makes  wires  along  with  contact 
bonds,  a  significant  amount  of  gold  may  be  used  between  all  of  the  samples.  If  the  gold  is 
activated  above  permitted  shipping  levels,  it  may  take  up  to  a  month  before  the  samples 
could  be  retrieved.  Indium  also  has  a  large  thermal  cross  section,  but  its  short  half-life 
means  that  its  activity  the  following  day  would  be  greatly  reduced.  However,  indium  has 
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a  melting  point  of  156°  C  (429  K),  which  makes  it  undesirable  for  armealing  to  200°  C.  A 
nickel  /  silver  wire  was  finally  chosen  for  wires  on  the  samples. 

Contact  Formation 

The  ability  to  make  effective  ohmic  contacts  for  silicon  carbide  is  a  major 
research  area.  All  of  the  benefits  of  a  high-temperature,  high-power  semiconductor  are 
lost  if  the  contacts  cannot  survive  the  environment.  Some  researchers  describe  the 
formation  of  good  contacts  as  critical  to  the  future  of  SiC  devices. 

Contacts  at  metal/SiC  interfaces  are  either  a  Schottky  (rectifying)  contact,  or  an 
ohmic  contact.  A  Schottky  contact  resembles  the  behavior  of  a  diode,  rectifying  the 
current  to  pass  preferentially  in  one  direction.  Ohmic  contacts  behave  like  resistors, 
increasing  current  flow  proportionally  with  increased  applied  voltage.  An  ohmic  contact 
is  produced  when  the  Schottky  Barrier  Height  (SBH)  is  minimized.  For  this  project,  it  is 
important  to  make  contacts  with  ohmic  behavior  and  low  resistance. 

For  a  long  time,  nickel  has  been  the  material  of  choice  for  good  ohmic  contacts  on 
n-type  silicon  carbide.  It  was  initially  believed  that  the  formation  of  Ni2Si  was  the  reason 
good  contacts  formed.  The  Ni  -  Si  reaction  separates  the  SiC  leaving  behind  unreacted 
carbon,  and  was  a  concern  for  use  in  a  high  temperature  environment  [31].  The  fear  was 
that  the  carbon  may  bond  with  another  material  and  degrade  the  contact.  Studies  show 
that  Ni2Si  forms  at  much  lower  temperatures  than  the  annealing  temperature  needed  for 
contact  formation  [29].  In  most  cases,  annealing  temperatures  above  900°C  are  needed  to 
make  low  resistance  ohmic  contacts,  with  the  exception  of  highly  doped  samples  [48]. 
One  study  notes  a  significant  improvement  in  resistance  by  armealing  to  1000°  C. 
However,  Ni2Si  forms  as  low  as  500°  C.  Studies  of  the  unreacted  carbon  suggest  that 
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graphitization  of  the  material  was  the  driving  foree  behind  ohmic  contact  formation  [30]. 
It  was  found  that  graphite  formed  at  temperatures  of  900°  C,  very  close  to  the  necessary 
annealing  temperatures.  It  was  also  found  that  materials  that  act  as  graphitic  catalysts, 
such  as  Ni  and  Co,  make  good  ohmic  contacts  on  n-type  SiC.  Lu  claims  that  Ni  is  used 
for  ohmic  contacts  on  n-type  material  not  because  we  understand  why  it  is  a  good  choice, 
but  only  that  it  works.  He  comments  that  Ni,  as  well  as  Co,  make  good  ohmic  contacts 
because  they  are  excellent  catalysts  for  graphite  formation  in  n-type  SiC  [30].  In  another 
paper,  Lu  points  out  the  silicides  were  once  believed  to  be  the  source  of  ohmic  contact 
formation,  but  the  annealing  temperature  to  form  a  silicide  does  not  coincide  with  the 
temperature  to  form  an  ohmic  contact  [29]. 

While  good  contacts  are  easily  fabricated  for  n-type  SiC,  it  is  very  difficult  to 
make  suitable  ohmic  contacts  for  p-type  SiC.  Successful  ohmic  contacts  have  only  been 
formed  on  highly  doped  (lO'*  cm‘^)  samples.  Ohmic  contacts  normally  form  by  bonding  a 
metal  possessing  a  work- function  equivalent  to  the  electron  affinity  of  the  semiconductor. 
However,  no  known  metal  possesses  a  work  function  equivalent  to  6eV,  the  electron 
affinity  for  p-type  SiC,  and  thus  a  Schottky  contact  will  always  form.  However,  contacts 
can  be  made  where  electrons  turmel  through  the  barrier  in  an  ohmic  manner.  For  the  past 
few  years,  much  research  has  been  conducted  on  making  effective  ohmic  contacts  for 
silicon  carbide. 

Finding  a  suitable  material  remains  a  difficult  task.  Aluminum,  a  common  dopant 
for  p-type  SiC,  was  recognized  as  having  electrical  properties  suitable  for  good  ohmic 
contacts.  Aluminum  has  a  low  melting  point,  which  is  not  suitable  for  use  in  high 
temperature  applications.  Titanium  was  added  to  aluminum  to  increase  the  melting  point. 
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The  AFIT  labs  have  traditionally  used  an  Al/Ti  alloy  eonsisting  of  90%  A1  and  10  Ti% 
(90/10)  for  SiC  contacts.  However,  a  study  on  optimizing  titanium  content  in  the  alloy 
found  that  using  a  70/30  mixture  improves  resistivity,  and  increases  reproducibility 
significantly  [10].  A  morphological  study  of  70/30  contacts  found  that  the  Al/Ti  alloy 
contained  a  liquid  phase  that  “spiked”  into  the  material  [39].  The  90/10  alloy  resulted  in  a 
few  long  spikes,  while  the  70/30  resulted  in  shorter  spikes  with  two  orders  of  magnitude 
more  abundance.  The  conclusion  of  the  study  was  that  the  addition  or  subtraction  of  a 
spike  in  the  90/10  alloy  lead  to  a  much  larger  variance,  and  explained  the  difficulty  of 
reproducing  results  with  that  mixture.  The  size  of  the  70/30  spikes  was  65  nm  maximum 
in  depth. 

Hall  Effect  Measurements 

In  1879,  as  a  master’s  student,  Edwin  Hall  discovered  the  phenomena  now  known 
as  the  Hall  effect.  Hall  noticed  that  when  a  conductor  was  exposed  to  a  magnetic  field 
oriented  perpendicular  to  the  current  through  the  conductor,  a  voltage  appeared  across  the 
height  of  the  conductor.  This  phenomenon  is  related  to  the  Lorentz  force,  as  all  of  the 
carriers  through  the  conductor  are  experiencing  a  qvxB  force  as  they  traversed  though 
the  device  as  depicted  by  Figure  2.  Hall  described  the  effect  as  “pressing  electricity.” 


14 


Figure  2:  Hall  effect  in  a  planar  bar.  A  current  is  run  perpendicular  to  an  external  magnetic  field, 
resulting  in  a  force  on  the  carrier.  Source  NIST  [41] 


While  Hall  used  metal  for  his  experiments,  Hall  effect  techniques  prove  very 
helpful  in  the  characterization  of  semiconductors,  and  shall  be  applied  in  this  experiment 
to  describing  the  mobility  of  neutron  irradiated  silicon  carbide. 

Initially,  the  Hall  effect  was  used  to  describe  the  resistance  of  a  particular 
material.  Eventually,  resistance  was  no  longer  considered  a  characteristic  of  a  particular 
material,  since  values  were  geometrically  dependent.  Resistance  was  replaced  as  a 
characteristic  property  by  resistivity.  Resistivity  quantifies  the  difficulty  of  a  carrier  to 
move  through  a  material  [41]. 

Resistivity  is  not  independent  of  other  material  properties,  and  varies  due  to 
changes  in  carrier  concentration  in  a  material.  Resistivity  (p)  is  now  denoted  as  a 
function  of  carrier  mobility  (p^),  carrier  density  (n),  and  carrier  charge  {q)  as  related  in 
equation  4. 

p  =  —  (4) 

fx^nq 
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Since  carrier  density  varies  with  temperature,  resistivity  is  also  implieitly  dependent  on 
temperature.  Thus,  for  this  study  it  is  important  to  charaeterize  temperature  dependent 
ehanges  in  resistivity  to  identify  damage  dependent  ehanges  in  resistivity. 

Hall  effect  measurements  now  aim  at  measuring  mobility  as  a  material’s  intrinsie 
electrieal  property.  For  eomparing  different  materials,  their  mobility  should  be  eompared 
[41].  Hall  effect  measurements  measure  the  Hall  mobility  of  the  material,  and  can  be 
related  to  the  conduetion  mobility  by  the  relation: 

“1-1%  (5) 

Mobility,  while  being  unique  to  a  material,  differs  based  on  the  crystal  plane  of  the 
eonduetion  medium.  Mobility,  thus  is  a  second  rank  tensor  differing  in  direetion  [51]. 
Since  mobility  is  only  eompared  as  a  sealar  value,  it  is  implied  that  this  value  is  an 
average  over  the  entire  crystal. 

There  are  two  common  methods  of  taking  Hall  effect  measurements,  the  Van  der 
Pauw  method,  and  the  Hall  bar  method.  The  Hall  bar  method  requires  a  special  geometry, 
while  the  Van  der  Pauw  method  can  be  used  with  almost  any  geometry.  Due  to  the  added 
eomplexity  of  creating  a  Hall  bar  geometry,  the  Van  der  Pauw  method  was  used  in  this 
study.  The  Van  der  Pauw  method  eombines  4  point  Van  der  Pauw  resistivity 
measurement  with  Hall  measurements  to  determine  mobility. 

Hall  measurements  are  particularly  suited  for  studying  surfaee  defeets.  By 
applying  a  strong  magnetie  field,  eharge  earriers  are  eonfined  to  a  thin  region  of  material. 
Since  pre-knowledge  of  the  depth  is  necessary,  a  thin  conduetion  layer  is  required,  and  is 
aecomplished  by  using  samples  possessing  a  thin  epitaxial  layer. 
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III.  Experiment 


This  study  was  accomplished  in  four  phases,  equipment  setup,  sample  fabrication, 
device  irradiation,  and  temperature  dependent  measurement. 


Equipment  Setup  and  Testing 


The  first  priority  of  this  project  was  setting  up  and  testing  the  measuring 
equipment.  The  Keithley  system  110  was  recently  moved  from  another  laboratory  in  the 
school.  While  minimal  documentation  was  included  with  the  system,  the  manuals  for 
most  components  have  been  located  on  their  manufacturers’  respective  website.  The 
system  components  are  listed  in  Table  3.  The  system  depicted  in  Figure  3,  contains  a 
Walker-Scientific  Electromagnet,  and  is  controlled  by  a  computer  program  over  a  general 
programming  interface  bus  (GPIB). 


Table  3:  Equipment  List  for  Keithley  System  110 


Make 

Model 

Name 

Keithley 

617 

Electrometer 

Keithley 

196 

Voltmeter 

Keithley 

— 

Buffer  Amplifier 

Keithley 

706 

Scanner 

Keithley 

220 

Current  Source 

Lakeshore 

DRC-91CA 

Coolant  Controller 

Walker-Scientific 

— 

Electromagnet 

Walker-Scientific 

MG-3D 

Gaussmeter 

Walker-Scientific 

— 

Current  Controller 

The  system  can  be  controlled  either  manually  or  by  computer  program.  Manual 
operation  requires  opening  and  closing  channels  on  the  scanner  in  order  for  a  particular 
device  to  have  access  to  the  sample  from  the  quad  buffer  amplifier.  Unfortunately, 
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documentation  for  which  channels  conneet  to  whieh  devices  exists  only  in  the  form  of  a 
circuit  diagram. 

A  computer  program  written  by  Keithley  in  QBasic  and  modified  at  AFIT  was 
supplied  with  the  system.  The  eomputer  eontroller  eommunieates  to  the  system  using  a 
general  programming  instrument  bus  (GPIB)  or  IEEE  488  bus. 


Figure  3:  Keithley  System  110 


The  measuring  equipment,  whieh  had  not  been  used  for  some  time,  required 
extensive  assembly  and  troubleshooting  to  return  to  working  order.  Due  to  the  time 
required,  troubleshooting  was  accomplished  in  parallel  with  fabrication  and  irradiation. 
The  greatest  impediment  was  the  lack  of  documentation  stating  ehanges  in  the  system 
from  its  factory  delivered  state.  Three  major  problems  are  worth  discussing  -  computer, 
buffer,  and  power  souree. 
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The  original  documentation  suggests  running  the  supplied  SYSTST  program 
when  first  setting  up  and  testing  the  system.  However,  all  programs  executed  with  the 
supplied  computer  failed  with  an  error  stating  the  bus  could  not  be  initialized.  This 
complicated  troubleshooting  because  the  system  manual  suggests  not  attempting  to 
manually  operate  the  system  until  familiarity  with  computer  controlled  operation  is 
achieved.  As  such,  the  manual  reads  as  a  programming  guide,  with  no  information  about 
how  to  go  about  manual  operation. 

Difficulty  launching  the  original  Keithley  program  lead  to  the  development  of  an 
in-house  Visual  Basic  based  program.  The  new  program  was  written  in  the  time  span  of  a 
month,  not  allowing  enough  time  for  error  checking  code  to  be  fully  implemented.  The 
code’s  modular  format  permits  future  users  to  easily  develop  new  interfaces  without 
having  to  redevelop  the  entire  package. 

The  original  design  of  the  System  110  supplied  power  to  the  quad  amplifier  buffer 
through  the  scanner.  The  buffer  had  been  modified  to  be  powered  with  an  external 
voltage  source.  The  required  power  source  must  be  able  to  supply  both  positive  and 
negative  polarity,  sourcing  -15V  low,  and  simultaneously  15V  high.  This  supplies  power 
to  the  unity  gain  amplifiers  required  for  high  resistance  testing.  Fortunately,  the  SiC 
devices  were  of  low  enough  resistance  to  use  the  low  resistance  terminals. 

After  successfully  measuring  a  voltage  drop  across  a  resistor  with  a  known 
applied  current,  the  magnet  was  to  be  tested.  Following  the  schematics  supplied  by  a 
Walker-Scientific  applications  engineer,  it  was  determined  that  a  coolant  flow  safety  was 
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not  electrically  closing  as  designed.  Permanently  shorting  the  circuit  until  the  proper  part 
could  be  ordered  offered  a  temporary  solution. 

Sample  Fabrication 

Fabrication  included  cutting,  cleaning,  and  forming  ohmic  contacts.  Cutting 
samples  into  an  optimal  geometry  can  minimize  errors  in  Van  der  Pauw  measurements. 
Cleaning  the  sample  removes  impurities  that  can  affect  the  creation  of  a  good  contact. 
While  good  contacts  can  be  made  easily  for  n-type  SiC,  it  is  very  difficult  to  make  good 
ohmic  contacts  for  p-type  SiC.  Wiring  was  connected  prior  to  irradiation  to  minimize 
heating.  Due  to  the  unique  circumstances,  of  this  experiment,  traditional  wiring 
techniques  were  not  acceptable,  and  proved  challenging. 

Sample  Cutting 

The  different  methods  for  Hall  effect  measurements  demand  different  geometries 
of  use.  The  Van  der  Pauw  method  used  in  this  experiment  can  be  executed  with  an 
arbitrary  geometry.  However,  in  order  to  use  derived  equations,  and  minimize  error,  there 
exist  recommended  geometries. 

NIST  recommends  use  of  a  “cloverleaf  ’  geometry  as  depicted  in  Figure  4. 
However,  lithography  methods  must  be  used  in  order  to  develop  this  complex  design. 
Satisfactory  results  can  be  obtained  using  a  square  geometry,  with  negligible  loss  in 
accuracy. 
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Figure  4:  Recommended  Geometries  from  NIST  for  use  with  Van  der  Pauw  style  Hall 

Measurements.  Source  [41] 


In  order  to  ensure  aceurate  data  a  sufficient  number  of  samples  must  be  measured. 
Jones’  results  from  a  previous  thesis  suggest  that  sample  deviation  for  n-type  samples 
was  only  3%  error  while  deviation  of  p-type  samples  was  about  10%  [20,  54-57].  Small 
sample  statistics  can  be  used  to  determine  the  minimum  number  of  samples  required  for 
this  study.  Equation  6  was  used  to  determine  the  minimum  number  of  samples  required 
for  the  desired  statistics. 


II 

where 

I 

confidence  interval. 

t 

t  student  distribution. 

V 

degrees  of  freedom  =  N- 1 , 

N 

number  of  samples,  and 

s 

sample  standard  deviation. 

(6) 


Since  the  n-type  SiC  has  a  small  variance,  only  4  samples  per  fluence  are  required  for 
95%  confidence  of  a  5%  interval.  P-type  SiC  has  a  much  larger  variance;  a  90% 
confidence  interval  with  a  resolution  of  10%  was  chosen  to  keep  the  sample  size  to  a 
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manageable  5  samples  per  fluence.  Over  the  range  of  4  fluences,  20  samples  from  each 
wafer  will  satisfy  requirements  for  this  thesis. 

Several  wafers  of  4H-SiC  were  purchased  from  Cree  Inc  in  support  of  a  previous 
thesis.  Enough  previously  cut,  but  otherwise  untouched,  remained  from  that  experiment 
to  be  used  with  this  study  with  satisfactory  statistics.  Figure  5  shows  a  5mm  x  5mm 
sample  with  a  penny  for  scale. 


Figure  5:  Example  5mm  x  5mm  Sample. 

Note  that  sample  is  roughly  the  size  of  Lincoln's  head  on  a  penny. 

A  few  samples  had  to  be  cut  by  a  wafer  saw  to  get  the  desired  5mm  x  5mm  square 
dimensions.  On  the  day  of  the  cut,  the  wafer  saw  did  not  completely  cut  through  the 
material,  producing  only  a  deep  grove.  Excess  material  was  removed  by  being  snapped  at 
the  groove.  While  this  procedure  is  common  for  silicon,  it  can  lead  to  cracking  in  silicon 
carbide  due  to  a  non-planar  crystal  boundary.  Figure  6  depicts  the  groove  left  on  one 
sample,  and  a  crack  caused  on  another  sample  due  to  snapping  off  the  extra  material. 
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Figure  6:  (Left)  A  groove  cut  into  a  sample.  (Right)  A  crack  caused  by  snapping  off  excess  material. 


While  enough  pre-cut  samples  existed  to  meet  minimum  desired  statistics,  a 
certain  amount  of  attrition  should  have  also  been  factored  into  the  statistics.  Over  the 
course  of  the  project,  several  samples  were  either  lost,  did  not  display  ohmic  behavior,  or 
were  damaged  before  a  measurement  could  be  taken.  As  a  result,  the  number  of  samples 
measured  was  less  than  the  number  required  as  calculated  above.  After  ensuring  that  the 
samples  possessed  the  desired  geometry,  the  samples  were  cleaned  in  preparation  for 
deposition  of  the  ohmic  contacts. 


Device  Cleaning 

The  purpose  of  cleaning  is  to  remove  any  impurities  that  may  have  occurred  from 
handling  or  chemical  interaction  from  the  atmosphere  in  order  to  create  a  good  bond 
between  the  metal  contact  and  the  semiconductor.  Cleaning  methods  vary  depending  on 
purpose  and  preference,  however  common  elements  are  present. 

The  selected  cleaning  process  consisted  of  a  rinse  of  TCE,  acetone,  methanol,  and 
deionized  (DI)  water,  followed  by  N2  Dry.  The  sample  was  held  at  the  edges  by  tweezers 
instead  of  being  held  by  a  spinner  due  to  the  small  sample  size.  Rinsing  took  a  period  of 
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10  to  12  seconds  per  solvent  by  counting  slowly  to  12.  Care  had  to  be  taken  when  drying 
with  the  N2  gun  since  the  gas  velocity  was  high  enough  to  blow  the  sample  across  the 
hood.  One  sample  was  lost  in  this  fashion.  The  samples  were  then  diamond  scribed  with 
their  respective  positions  on  the  original  wafer,  and  cleaned  again.  Figure  7  shows  the 
chemicals  used  in  the  chemical  hood. 


Figure  7:  Cleaning  Supplies  in  Hood.  Samples  were  eleaned  prior  to  eontact  deposition 

Additionally,  the  quartz  tube  used  in  the  tube  furnace  also  required  cleaning. 

Since  the  material  previously  used  in  that  tube  was  unknown,  a  strong  etch  had  to  be 
performed  using  the  acid  aqua  regia.  Aqua  regia,  a  3:1  mixture  of  HCLiHNOs,  is  a  strong 
oxidizing  agent  which  gets  its  name  (royal  water)  from  the  ability  to  dissolve  the  royal 
metals  such  as  gold  and  platinum.  The  etching  process  ensures  that  residue  left  on  the 
tube  from  previous  anneals  of  different  materials  does  not  contaminate  current  samples. 
The  tube  was  left  filled  with  acid  for  about  half  an  hour,  when  the  bubbling  from 
chemical  reaction  ceased.  The  tube  was  then  labeled  “SiC”  to  inform  future  users. 
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Forming  Ohmic  Contacts 


Ohmic  contacts  were  deposited  as  soon  as  all  samples  were  cleaned.  Contacts 
were  deposited  using  the  E-beam  system  depicted  in  Figure  8.  Contaets  were  placed  on 
the  periphery  of  the  samples,  as  close  to  the  comers  as  practical.  Layers  of  150  nm  were 
plaeed  on  the  samples.  Niekel  was  used  for  n-type  samples,  and  A1  /Ti  (70/30)  was  used 
for  p-type  material.  Since  similar  Hall  effeet  studies  are  being  done  at  AFIT,  a  suitable 
contact  mask  already  existed. 


Figure  8:  E-Beam  Deposition  System.  This  system  deposited  the  metal  contacts  onto  the  samples 
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As  an  initial  check  of  contact  deposition,  samples  were  examined  by  microscope 
inside  of  the  clean  room.  Figure  9  is  an  example  of  a  sample  image  captured  by  computer 
from  the  microscope.  The  samples  were  so  large  that  at  the  microscope’s  broadest  setting, 
the  entire  sample  could  not  be  captured.  The  image  collected  from  each  sample 
corresponds  to  the  worst  looking  comer. 


Figure  9:  Sample  under  microscope 


After  all  contacts  were  deposited  on  all  samples,  they  were  annealed  in  an 
annealing  furnace.  Before  samples  could  be  annealed,  the  quartz  tube  had  to  be  put 
through  a  heating  cycle  to  ensure  that  any  residual  moisture  was  released  from  the  tube 
since  it  had  just  been  cleaned. 

The  annealing  procedure  was  as  follows.  First,  the  sample  was  loaded  into  the 
tube,  but  outside  of  the  furnace.  Figure  1 1  depicts  samples  in  the  staging  area  outside  of 
the  furnace.  Next,  a  supply  of  argon  gas  was  connected  to  the  tube.  While  nitrogen  was 
originally  selected,  argon  was  used  during  actual  aimealing  because  nitrogen  is  a  known 
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dopant  of  SiC,  and  argon  is  also  an  inert  gas.  This  had  the  advantage  of  inert  gas  flow 
without  aoeidentally  ehanging  the  doping  level. 


Figure  10:  Samples  annealing  in  furnace.  Samples  were  annealed  at  temperatures  of 
900°  C  andll00°  C  for  p-type  and  n-type  respectively. 

Once  inert  gas  was  flowing,  the  furnace  was  turned  on.  It  usually  took  about  20 
minutes  for  the  furnace  to  reach  the  set  temperature.  Once  the  required  temperature  was 
reached,  the  samples  were  inserted  into  the  hot  furnace  for  the  prescribed  period,  and  then 
removed  to  the  staging  area.  N-type  samples  required  an  annealing  time  of  5  minutes  at 
1100°  C,  and  p-type  samples  required  2  minutes  at  900°  C.  To  prevent  oxidation,  the 
samples  were  allowed  to  remain  under  argon  gas  flow  while  the  samples  and  the  furnace 
cooled  down. 
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Figure  11  Samples  cooling  after  anneal 

The  effectiveness  of  the  annealing  process  was  measured  using  a  probe  station  to 
take  IV  data  of  each  sample.  Figure  12  shows  an  image  of  the  probe  station  used  to  test 
the  newly  formed  ohmic  contacts.  A  voltage  was  applied  to  contact  1  &  3,  measuring  the 
resultant  current.  The  process  was  repeated  with  contact  2  &  4  on  the  sample. 


Figure  12:  Semi-Conductor  Characterization  Probe  Station.  This  station  was  used  to  take  IV 

measurements. 
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A  properly  fabricated  contact  is  one  with  low  resistance  ohmic  behavior.  Most 
samples  displayed  the  proper  ohmic  behavior  with  linear  IV  plots.  An  example  of  the 
resulting  IV  curve  from  a  good  sample  is  Figure  13.  Figure  14  depicts  the  result  from 
contacts  exhibiting  non-ohmic  behavior. 


Voltage  [V] 

o  C9N  o  C9NN -———-Linear (C9NN) - Linear (C9N) 

Figure  13:  IV  curve  for  n-type  sample  C9.  Note  the  linear  fit. 
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I-V  Curve  for  Sample  H3N 


Figure  14:  IV  curve  for  n-type  sample  H3.  Note  the  rectifying  behavior. 


Wire  Bonding 

Two  methods  are  common  for  connecting  square  Van  der  Pauw  samples  to  a  Hall 
system.  The  first  method  involves  using  probes  to  press  against  the  ohmic  contacts.  The 
second  method  is  to  connect  thin  wires  with  a  solder,  usually  indium.  The  Hall  effect 
system  used  in  this  research  uses  the  latter  method. 

Typically,  the  Hall  wires  are  soldered  to  the  sample  holder.  This  experiment 
hinges  on  minimizing  premature  heating  of  the  sample,  so  post  irradiation  soldering  is 
avoided  as  much  as  possible.  To  avoid  soldering  wires  post  irradiation  (and  inadvertently 
heating  the  sample),  wires  were  placed  on  the  samples  prior  to  irradiation.  Figure  15  is  a 
picture  of  a  sample  with  attached  wires. 
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Figure  15:  Sample  with  attached  wires 

Since  pre-irradiation  attaehment  places  additional  materials  in  the  reactor, 
connecting  wires  are  required  to  posses  a  low,  fast  spectrum  neutron  capture  eross 
sections,  and  survive  environments  of  temperature  ranging  from  77  K  to  500  K. 

Common  materials  used  in  electronies  inelude  gold,  indium,  silver,  aluminum, 
and  copper.  Copper  and  gold  were  avoided  due  to  activation.  Soldering,  wire  bonding, 
and  conduetive  epoxy  were  the  three  bonding  methods  available. 

I  initially  settled  on  wire  bonding  aluminum  wires  to  the  samples.  The  thin  wires 
used  for  wire  bonding  would  have  had  the  advantage  of  low  mass,  and  thus  low 
aetivation.  Furthermore,  aluminum  is  a  low-aetivation  material,  with  aeeeptable  electrieal 
properties.  Thus  I  chose  to  wire  bond  aluminum  wires  to  the  ohmie  pad.  This,  however 
presented  additional  problems  with  handling  since,  bonded  wires  are  not  usually  intended 
to  be  handled  in  the  manner  proposed,  and  risked  snapping  off. 

Soldering  was  then  explored.  The  diffieulty  with  using  solder  is  its  generally  low 
melting  point.  Several  high  temperature  lead  /  indium  solders  were  located  at  AFRL. 
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Unfortunately,  none  of  the  solders  properly  whet  the  material,  not  sticking  properly  to 
either  the  nickel  or  AlTi  contacts.  Conductive  epoxy  was  the  last  remaining  option. 

The  conductive  epoxy  is  a  liquid  containing  silver  flakes  that  highly  concentrate 
when  cured.  A  silver  /  nickel  wire  was  selected  to  compliment  the  epoxy. 

Since  the  wires  would  be  irradiated,  material  selection  was  important  due  to 
activation  concerns.  To  minimize  activation,  the  wires  were  cut  to  just  long  enough  to 
reach  the  posts  on  the  sample  holder.  Initial  attempts  to  mount  samples  on  the  holder 
showed  that  the  epoxy  connection  was  fragile.  The  wires  were  then  bent  at  right  angles  to 
relieve  stresses  found  at  the  wire  /  contact  junction.  Figure  16  shows  a  sample  with 
shortened  bent  wires. 


Figure  16:  Sample  with  bent  wires  in  preparation  for  irradiation. 

While  this  configuration  was  more  resistant  to  wire  detachment,  failures  still 
occurred. 

Sample  Irradiation 

Samples  of  4H-SiC  were  exposed  to  fiuences  of  fast  neutrons  ranging  from  lO'"^ 

TV  *  CTTl 

to  10*^ - ^  at  the  Ohio  State  University’s  Nuclear  Research  Reactor  (OSUNRR). 

cm 

Corresponding  irradiation  times  are  calculated  in  Appendix  B.  Due  to  the  high  fiuence 
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and  required  dose,  the  Central  Irradiation  Faeility  (CIF)  was  used  for  irradiation.  The  CIF 
is  a  tube  that  allows  aceess  to  the  center  of  the  reactor.  The  reactor’s  rabbit  tube  was  also 
an  option  for  neutron  irradiation,  but  was  not  used  since  the  pneumatic  tube  has  a 
tendency  to  shock  its  contents  with  sudden  jolts  that  may  have  damaged  the  wire 
connection. 

In  order  to  ensure  low  activation,  the  samples  were  cadmium  shielded.  Cadmium 
is  a  unique  material,  allowing  most  of  the  fast  neutron  spectrum  to  pass  through,  while 
filtering  most  of  the  thermal  neutrons. 


012345678 

Neutron  Energy  [eV] 


Figure  17:  Cadmium  neutron  absorption  cross  section  as  a  function  of  energy.  Cadmium  has  the 
unique  property  of  a  large  thermal  cross  section  and  a  small  fast  crossection  [24] 


Since  thermal  neutrons  do  not  have  the  energy  to  create  the  damage  desired  in  this 
experiment,  using  cadmium  shielding  will  not  skew  results  of  irradiation. 
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Figure  18:  Sample  and  its  Container  prior  to  irradiation. 


The  SiC  devices  were  placed  in  plastic  containers  just  big  enough  to  fit  the 
samples.  Figure  18  shows  the  container  for  a  device,  while  Figure  19  shows  a  stack  of 
several  samples  in  a  run  before  applying  the  cadmium  wrapping.  It  should  be  noted  that 
the  amount  of  cadmium  required  for  6  samples  exceeded  allowable  amounts  of  negative 
reactivity  for  the  reactor.  To  prevent  reactivity  problems,  only  3  samples  were  irradiated 
at  a  time,  reducing  the  amount  of  Cd  shielding.  For  higher  doses,  two  samples  were 
placed  in  one  container  to  further  reduce  the  amount  of  cadmium  required  to  cover  all 
samples  of  the  run.  The  required  thickness  of  cadmium  can  be  calculated  with  (7). 

llk=e-^‘  (7) 

,  ^  -InlO.lOl  2.250 

leads  to,  t  = - ^ ^  « - 

2  2 


where. 

llh  ~ 

intensity. 

2 

macroscopic  crossection,  and 

t 

thickness. 

Using  2=463  (1/cm),  and  90%  thermal  neutrons  filtered,  t  =  0.005  cm  or  2  mil. 
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While  the  eadmium  wrapping  prevented  activation,  a  disadvantage  of  the 
cadmium  became  apparent.  The  cadmium  acted  as  a  poison,  a  neutron  absorbing 
substance  adding  negative  reactivity.  To  counteract  the  poison,  control  rods  were 
withdrawn.  For  safety  reasons,  control  rod  withdrawal  is  limited.  Due  to  additional 
negative  reactivity,  the  number  of  samples  stacked  were  reduced  from  6  to  3  in  order  to 
reduce  required  cadmium  shielding.  For  future  experiments,  it  is  suggested  to  combine 
samples  together  in  a  small  space  to  reduce  the  cadmium  requirement.  Additionally, 
using  the  minimum  thickness  of  cadmium  required  would  reduce  negative  reactivity. 
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Figure  20:  Containers  wrapped  in  cadmium.  Cadmium's  high  thermal  neutron  cross  section  helps  to 

minimize  activation  of  the  samples. 

Measurement  Procedure 

Once  the  system  was  operational,  a  method  had  to  be  devised  to  mount  the 
samples  to  the  eold  head.  The  nature  of  this  experiment  requires  the  sample  to  stay  as 
cold  as  possible  during  mounting,  eliminating  use  of  solders,  which  would  add  heat 
through  the  wires.  Additionally,  the  junetion  between  the  wire  and  the  ohmic  contact, 
which  was  held  by  silver  epoxy  proved  to  be  very  fragile. 

Several  mounting  techniques  were  attempted  without  satisfaetory  results.  These 
techniques  ineluded  using  conductive  epoxies,  solders,  and  attaching  thinner  wires.  The 
method  finally  implemented,  required  removing  the  pre-conneeted  wires,  and  connecting 
a  new  wire  (whieh  is  soldered  onto  the  sample  mount)  to  the  sample  using  a  conductive 
silver  epoxy.  Since  euring  the  epoxy  was  not  an  option  after  irradiation,  a  dab  of  rubber 
cement  was  plaeed  on  the  wire  to  hold  it  in  place.  Using  an  uneured  epoxy  added 
minimal  resistance  to  the  contacts. 
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After  ohmic  contacts  had  been  applied  to  the  surface  of  the  SiC,  the  initial 
mobility  of  the  samples  needed  to  be  recorded.  The  same  process  was  followed  for  the 
Hall  effect  measurements  before  and  after  irradiation. 

Samples  were  characterized  after  device  fabrication  and  irradiation.  Three 
samples  per  type  were  reserved  for  initial  characterization,  corresponding  to  a  neutron 
fluence  of  zero.  The  remaining  12  samples  were  irradiated  at  OSUNRR  facility  as 
described  in  the  irradiation  section. 

In  order  to  characterize  from  lOOK  to  340K,  the  unirradiated  samples  must  first  be 
cooled  with  liquid  nitrogen.  Irradiated  samples  are  cooled  in  liquid  nitrogen  to  damage 
since  vacancy  and  interstitial  damage  are  immobile  at  temperatures  less  than  80  K 
[17:82]. 

Samples  were  bonded  to  the  sample  mount  using  rubber  cement.  The  wires  on  the 
sample  were  then  soldered  onto  the  pins  of  the  sample  mount.  Finally,  the  sample  mount 
was  replaced  into  the  cold  head. 

Hall  measurements  over  a  temperature  range  from  lOOK  to  340K  at  lOK 
increments  were  measured.  The  lower  end  was  decided  due  to  limitations  on  keeping  the 
device  at  liquid  nitrogen  temperature  during  mounting,  accepting  a  23K  buffer  above  77 
K  as  reasonable.  The  upper  limit  temperature  was  limited  by  the  bounds  of  the  Hall  effect 
system.  Measurements  reaching  temperatures  of  500K  were  desired,  but  a  system  capable 
of  both  low  and  high  temperature  measurements  was  not  available. 

Hall  Effect  measurements  are  widely  used,  and  Standard  F76  is  published  by 
ASTM  to  ensure  continuity.  However,  a  copy  of  the  standard  was  not  purchased  for  this 
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study.  NIST  has  a  recommended  procedure  published  on  there  website,  and  Runyan 
makes  additional  recommendations,  both  based  on  ASTM  F76  Standard  [41,  51].  A 
detailed  description  of  how  to  take  Hall  Effect  measurements  can  be  found  in 
Appendix  A 
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IV.  Results 


The  results  of  this  experiment  are  broken  down  into  two  seetions.  First,  the  SI  10 
was  unable  to  match  parameter  values  of  the  unirradiated  characterization  samples  with 
AFRL’s  Hall  system.  Characterization  data  is  compared  and  differences  between  systems 
is  identified  to  narrow  down  the  problem,  but  no  particular  feature  is  identified  as  the 
cause.  Second,  results  pertaining  to  irradiated  samples  are  reported.  No  features  seen  in 
the  data  suggest  annealing  behavior  below  a  temperature  of  3  5  OK.  Characteristics 
expected  from  neutron  damage  are  identified.  An  anomalous  spike  in  values  is  seen  on 
occasion  at  temperatures  of  about  280K,  and  explained  as  a  result  of  using  an  uncured 
epoxy  on  some  samples. 

The  following  data  was  collected  from  n-type  samples  irradiated  with  lO'^  and 

lO'^ - 1  MeV  equivalent  fluences,  along  with  an  unirradiated  set  of  samples.  P-type 

cm 

samples,  and  n-type  samples  irradiated  to  lO'"^  measured  due  to  time 

cm 

restrictions. 

Additionally,  the  temperature  range  of  interest  was  restricted  due  to  equipment 
limitations.  The  AFRL  Hall  system  had  a  ceiling  of  340K,  while  the  SI  10  system  had  a 
ceiling  of  425K.  Initially,  a  temperature  of  500  K  was  the  goal  of  the  experiment,  but 
measurements  could  only  be  recorded  up  to  each  system’s  respective  limit. 

A  few  samples  were  eliminated  due  to  the  inability  to  form  ohmic  contacts. 

Fifteen  out  of  twenty  contacts  of  each  type  developed  sufficient  contacts.  The  n-type 
samples  selected  had  contact  resistances  ranging  from  10  Q  to  40  Q,  which  lead  to 
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resistivities  on  the  order  of  10'^  Q-cm.  The  literature  states  that  good  ohmie  contaets 
should  have  resistivities  on  the  order  of  10'^  Q-cm.  P-type  samples  had  resistances  from 
90  Q.  to  200  Q,  but  were  not  measured. 

Sample  population  decreased  further  after  wire  contact  formation.  Due  to  the 
combination  of  stiff  wire  and  a  brittle  connection  to  the  ohmic  contact,  the  silver  epoxy 
often  separated  from  the  contact.  Samples  broken  prior  to  irradiation  were  removed  from 
the  pool  of  samples. 

SllO  &  Characterization  Measurements 

The  first  objective  of  this  experiment  was  to  set  up  and  characterize  a  Keithley 
SllO  Hall  effect  system.  Temperature  dependent  Hall  measurements  of  unirradiated 
samples  were  taken  using  AFRL/SN’s  Hall  system  to  compare  to  the  results  from  the 
SllO  system. 
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Figure  21:  Comparison  of  Resistivity  Measurements  between  SllO  Hall  System  and  AFRL’s  Hall 

System 

Measurements  of  resistivity  from  the  SllO  were  mueh  greater  than  those  taken  by 
AFRL  as  depleted  by  Figure  21.  This  result  is  surprising  beeause  the  most  signifioant 
differenee  between  the  two  systems  was  with  the  eonnection  of  the  sample  wires  whieh 
ean be  aceounted  for  in  the  ealculations.  The  SllO  samples  had  a  eured  silver  epoxy 
eonneetion,  while  an  uneured  silver  epoxy  was  used  on  the  AFRL  system.  Since  the 
cured  sample  has  an  increased  density  of  the  silver  flake  in  the  epoxy,  it  is  reasonable  to 
expect  the  cured  sample  to  result  in  a  lower  resistance.  Studies  with  irradiated  samples, 
discussed  in  the  next  section,  composed  of  both  cured  and  uncured  contacts  show  no 
difference  in  resistivity  results. 
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Figure  22:  Comparison  of  Hall  Mobility  between  SllO  and  AFRL  Hall  systems 

Figure  22  presents  the  results  of  Hall  mobility  for  both  systems,  and  is  ealeulated 
from  the  results  of  the  Van  der  Pauw  resistivity  and  earrier  eoneentration  measurements. 
In  general,  mobility  inereases  with  temperature  at  similar  rates  for  both  systems.  The 
figure  indieates  a  higher  mobility  for  the  SllO  measurements  than  those  measured  at 
AFRL.  Mobility  values  are  eharaeterized  on  the  order  of  10^  [om^/(Vs)]  using  the  AFRL 
system. 

Figure  23  depiets  the  mobility  results  from  the  SllO  system  alone.  The  literature 
reports  a  maximum  in  mobility  just  below  room  temperature.  Figure  23  shows  a  loeal 
maximum  at  280K  whieh  is  diffieult  to  diseem  due  to  the  high  varianee  in  the  mobility 
values.  This  peak  mobility  has  an  average  value  of  826  om^/(Vs).  This  peak  is  not 
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obvious  from  the  results  taken  at  AFRL.  A  seeond  order  polynomial  fit  is  shown  as  a 
rough  approximation  to  emphasize  the  general  trend.  It  is  not  a  good  fit  sinee  the 
maximum  is  a  result  of  a  ehangeover  in  driver  from  impurity  seattering  to  aeoustic 
phonon  deformation  scattering  [22], 
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Figure  23:  Characterization  of  Hall  Mobility  using  SllO  Hall  Effect  System 

Bulk  carrier  concentration  as  measured  by  the  two  systems  also  greatly  varies. 
Figure  24  suggests  a  linear  increase  in  carrier  concentration  with  temperature  as 
measured  by  the  SllO  system.  Theory  models  the  temperature  dependence  of  carrier 
concentration  as  an  exponential. 
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Figure  24:  Bulk  Carrier  Concentration  of  4H-SiC  using  SllO  Hall  System 


Figure  25  shows  much  greater  bulk  carrier  concentration  from  the  AFRL  system 
than  with  the  SllO  system.  The  difference  in  values  may  be  a  result  of  the  difference  in 
magnetic  fields  used  between  the  two  systems.  The  variance  in  carrier  concentration 
carries  over  into  the  mobility  calculation  and  increases  with  temperature  due  to  thermal 
noise.  Increasing  the  current  flow,  discussed  in  the  next  section,  dramatically  decreased 
the  amount  of  variance  by  adding  carriers  thus  increasing  Hall  voltage. 
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Figure  25:  Comparison  of  Bulk  carrier  concentration  between  SI  10  and  AFRL  systems 


As  a  tool  to  troubleshoot  the  program  for  SI  10,  a  visual  basic  program  was 
written  to  take  measured  currents  and  voltages,  and  calculate  resistivity,  mobility,  and 
carrier  density  using  the  same  routines  that  the  SI  lO’s  program  used.  Using  the 
measurements  taken  at  AFRL,  I  confirmed  that  the  program  correctly  calculated  the 
desired  parameters.  Figure  26  is  a  depiction  of  the  program  confirming  AFRL’s  results 
from  the  measurements. 

By  ruling  out  an  error  in  calculation,  the  source  of  the  discrepancy  can  be 
narrowed  down  to  either  an  error  in  taking  a  measurement,  or  a  difference  in  a  system 
parameter,  such  as  magnetic  field  strength.  In  the  process  of  setting  up  the  equipment, 
each  instrument  was  checked  many  times  against  a  known  load.  Additionally,  delays 
were  programmed  into  the  code  to  allow  a  device  to  settle  before  requesting  data  to 
ensure  proper  measurement.  Therefore,  I  suspect  that  the  difference  in  results  comes  from 
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a  difference  in  system  operation.  Differences  were  noted  in  the  quality  of  the  wire 
connection  to  the  sample,  the  strength  of  the  magnetic  field,  and  the  medium  filling  the 
cold  head. 


Figure  26:  Interface  for  Hall  effect  calculation  program. 


For  the  SI  10,  each  sample  had  the  original  wire  with  the  cured  epoxy  per  the 
original  plan.  Through  mounting,  unmounting,  and  handling,  many  of  the  fragile  wire 
connections  broke.  At  AFRL,  an  uncured  silver  epoxy  was  applied  to  the  exposed  ohmic 
contact  to  make  the  connection.  The  resistance  of  the  cured  samples,  as  measured  on  a 
digital  multimeter,  was  ~10  Ohms  while  the  uncured  samples  were  -50-100  ohms.  While 
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the  uncured  epoxy  had  a  slightly  higher  resistance,  it  wasn’t  enough  to  cause  the 
observed  differences. 

The  difference  in  magnetic  field  is  also  unlikely  to  have  played  a  significant  role. 

I  was  unable  to  get  the  SI  lO’s  magnet  above  IkG,  and  operated  it  at  0.7  kG.  AFRL 
operated  at  10  kG.  Since  the  magnetic  field  operates  during  the  Hall  measurement,  and 
not  the  Van  der  Pauw  measurement,  it  would  only  effect  the  measured  Hall  voltage,  and 
not  the  resistivity.  Likewise,  the  calculation  of  carrier  concentration  from  the  Hall  voltage 
already  includes  the  strength  of  the  magnetic  field,  so  even  with  different  fields,  the 
average  value  should  not  differ  as  much  as  observed  between  the  two  systems.  A  stronger 
magnetic  field,  however,  more  effectively  separates  the  charge  carriers,  increasing  Hall 
voltage,  and  decreasing  variance. 

The  SI  10  system  evacuated  the  cold  head  before  cooling  to  prevent  problems 
from  moisture  freezing.  The  AFRL  system  employs  a  stream  of  liquid  hydrogen  to 
displace  the  air  and  have  a  similar  effect.  The  advantage  of  a  streaming  gas  is  the 
reduction  in  cooling  time,  thus  increasing  the  sampling  rate.  It  is  unlikely  that  this 
difference  had  any  effect  on  the  system  as  well,  but  is  mentioned  for  completeness. 

Unfortunately,  this  qualitative  discussion  does  not  quantitatively  address  the 
accuracy  of  the  SI  10  system.  To  better  troubleshoot  the  SI  10  Hall  effect  system,  a 
standard  sample  should  be  used.  In  place  of  a  calibrated  standard,  a  device  similar  to  the 
ones  fabricated  in  this  thesis  made  of  silicon  could  be  used,  since  the  electrical  properties 
of  silicon  are  better  understood  than  SiC. 
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Radiation  Effects 

In  order  to  assess  the  effects  of  neutron  damage,  it  was  required  that  samples  be 
mounted  quickly  to  minimize  warming.  The  SI  10  did  not  allow  for  quick  sample  cooling, 
so  measurement  of  the  irradiated  samples  also  took  place  at  the  AFRL  sensors 
directorate.  The  SI  10  took  3  hours  to  cool  the  sample  to  100  K,  while  the  AFRL  system 
had  the  sample  back  to  100  K  in  30  minutes,  in  addition  to  the  mounting  time.  The 
characteristics  of  concern  are  resistivity.  Hall  mobility,  and  carrier  concentration.  Since 
all  of  the  irradiated  samples  were  measured  using  AFRL’s  equipment,  the  AFRL 
characterization  data  discussed  in  the  previous  section  is  used  for  comparison. 

Irradiated  samples  show  increased  resistivity  from  the  unirradiated  samples.  The 
average  results,  depicted  in  Figure  27,  suggest  that  resistivity  for  the  two  sets  of  irradiated 
samples  differ  from  the  characterization  set  by  a  multiplication  factor  in  agreement  with 
theory.  Theory  states  that  the  resistivity  of  an  irradiated  sample  will  be  the  resistivity  of 
the  unirradiated  sample  times  a  quantity  determined  by  the  neutron  fluence  received 
[37:242].  Resistivity  in  all  samples  appear  to  decrease  as  T"'  However,  a  difference 
between  the  resistivity  of  the  two  irradiated  sets  is  not  obvious  due  to  variance.  The 
lowest  resistivity  achieved  was  5.01E-5  Q-cm  occurring  at  3 3 OK,  supporting  the  claim 
that  good  n-type  ohmic  contacts  formed. 
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Figure  27:  Average  Resistivites  of  n-type  samples  by  neutron  fluence  received 

Up  to  340K,  the  decrease  in  resistivity  with  temperature  was  expected.  Above  this 
temperature,  resistivity  is  expected  to  increase  as  phonon  scattering  becomes  dominant 
over  Coulombic  scattering  from  impurities  and  thermal  excitation  of  carriers  into  the 
conduction  band. 

Note  that  an  anomalous  point  occurs  at  278K,  just  above  the  melting  point  of 
water  at  273K.  This  effect  was  only  observed  with  samples  measured  at  AFRL,  so  it  is 
not  likely  particular  to  the  material.  Also,  it  was  observed  with  irradiated  and  non- 
irradiated  samples  alike,  so  it  is  not  an  artifact  of  irradiation.  A  likely  explanation  is  that 
water  is  a  component  of  the  paste  being  used,  and  changed  the  resistance  of  the  contact  as 
it  melted  and  settled. 

It  is  important  to  note  that  sample  D6N  was  not  included  in  the  average  of 
irradiated  samples.  Sample  D6N  was  irradiated  with  a  fluence  of  lO'^  1-MeV  equivalent 
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neutrons,  but  did  not  show  any  signs  of  damage.  The  resistivity  of  D6N  in  Figure  28  is 
very  similar  to  the  average  resistivity  of  the  non- irradiated  samples  shown  in  Figure  27. 
For  this  reason,  D6N  was  not  included  in  the  average  of  the  lO'^  fluence  group. 


Temperature  [K] 


o  D6N  □  E7N  A  F4N 

Figure  28:  Resistivities  of  4H-SiC  samples  irradiated  with  10**  IMeV  equivalent  neutrons 


Carrier  density  decreases  with  irradiation  due  to  the  creation  of  active  acceptor 
defects  within  the  band  gap.  These  defects  rob  the  conduction  band  of  carriers, 
decreasing  carrier  density.  The  defect  is  formed  when  dislocated  silicon  or  carbon  atoms 
chemically  bond  with  an  impurity,  such  as  a  dopant,  within  the  lattice  [37:197].  Figure  29 
shows  almost  an  order  of  magnitude  decrease  in  carrier  density  to  irradiated  samples. 
Once  again,  the  difference  in  effect  between  the  two  fluences  is  small.  Carrier 

1  1  n^r.20  carriers .  .  ,  , 

concentrations  were  on  the  order  ot  10  - ^ — tor  irradiated  samples. 
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Figure  29:  Average  Bulk  Carrier  Density  for  n-type  4H-SiC  by  neutron  irradiation 


The  order  of  magnitude  (10^*’- 10^^  em‘^)  of  the  bulk  earrier  density  compared  to 
the  dopant  concentration  of  cm‘^  suggests  that  the  measured  value  is  high.  The 
likely  source  of  error  is  rooted  in  the  assumption  that  current  is  isolated  to  the  upper  most 
5  pm  epilayer.  The  n-type  samples  contained  an  n-type  buffer  region  lO'*  in 
concentration  and  0.5  pm  in  depth.  If  the  current  is  penetrating  further  into  the  substrate 
layer,  then  the  extra  depth  of  400  pm  can  significantly  alter  the  carrier  concentration 
calculation.  The  doping  level  of  the  substrate  is  unknown,  but  estimated  to  be  ~10'^. 
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Figure  30:  Average  Hall  Mobility  of  n-type  4H-SiC  separated  by  neutron  fluence 

A  decrease  in  mobility  with  irradiation  was  expected.  However,  Figure  30  shows 
that  the  difference  in  mobilities  between  irradiated  and  non-irradiated  sets  are  too  small 
to  be  quantified,  since  most  of  the  samples  lie  within  a  standard  deviation  of  each  other. 
This  suggests  that  neutron  damage  effects  carrier  concentration  more  than  mobility. 

The  largest  difficulty  in  analyzing  the  data  from  this  experiment  has  been  due  to 
the  large  variance  found  in  carrier  concentration.  It  was  found  that  increasing  the  sourced 
current  reduces  variance.  Sample  E9N  was  run  first  with  a  1mA  source,  and  again  with  a 
10mA  source.  The  results  from  the  second  case  show  much  reduced  variance  up  to  a 
temperature  of  200K,  where  thermal  noise  becomes  apparent.  This  improvement  is 
critical  to  identifying  features  in  the  data. 
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Figure  31  shows  that  the  inereased  eurrent  had  little  effeet  on  resistivity.  The 
benefit  is  strongly  eoupled  to  measurements  and  ealeulations  dependent  on  Hall  voltages. 
Figure  32  shows  features  in  earrier  eoneentration  that  were  not  previously  visible  due  to 
large  varianee.  A  minimum  earrier  eoneentration  is  identified  at  130K. 


Temperature  [K] 

o  1  mA  -  10  mA 


Figure  32:  Carrier  concentration  of  sample  E9N.  The  first  run  used  a  1mA  source  current.  The 

second  run  used  10mA. 

Due  to  the  improved  variance  on  carrier  concentration,  the  variance  on  mobility 
also  decreases.  Figure  33  shows  features  in  mobility  not  visible  in  previous  runs.  A 
maximum  mobility  appears  at  290K,  similar  to  the  supposed  maximum  described  in  the 
SI  10  section. 
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Since  the  inereased  current  test  was  only  performed  on  one  sample,  the  variance 
cannot  be  caleulated.  It  is  apparent  however,  that  choosing  a  proper  source  eurrent  is  very 
important  to  aceurately  measuring  a  Hall  voltage. 
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V.  Conclusions 


Concluding  Remarks 

This  study  attempted  to  identify  low-temperature  annealing  of  neutron  irradiated 
4H-SiC.  No  features  suggesting  annealing  below  340K  were  witnessed  in  this  study.  It 
was  determined  that  a  proper  souree  current  is  vital  to  identifying  features  that  may 
become  hidden  in  the  variance.  Since  studies  conflict  on  whether  damage  occurs  at 
temperatures  below  473K,  there  still  may  be  annealing  points  between  340K  and  500K. 
This  annealing  may  be  limited  to  the  surface  as  is  the  case  with  traditional  silicon. 

In  support  of  the  objective,  a  Hall  effect  system  was  set  up.  The  SI  10  system  set 
up  in  support  of  this  study  was  found  to  be  inadequate  due  to  a  long  cooling  time. 
Specific  problems  included  a  3  hour  cool  down  time,  providing  samples  with  too  much 
time  to  warm,  and  possibly  removing  damage  through  annealing.  Additionally,  a  probe 
based  sample  holder  would  be  beneficial  for  this  type  of  experiment  so  that  wire 
connections  could  be  avoided,  thus  reducing  activation  and  annealing  during  mounting. 

Silicon  carbide  samples  were  irradiated  and  characterized.  A  separate  set  of 
characterization  data  was  collected  from  both  the  in-house  SI  10  system  and  the  AFRL 
based  system.  The  scaling  of  each  result  differed  although  general  trends  appeared 
similar.  The  reason  for  the  different  values  may  lie  with  the  different  methods  of  wire 
connection,  in  combination  with  the  use  of  a  much  stronger  magnetic  field  at  AFRL.  The 
use  of  pastes  in  place  of  solders  is  not  recommended  since  some  components  may  freeze 
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and  melt  during  measurement.  Table  4  reviews  average  room  temperature  (300K) 
measurements  measured  using  both  systems. 


Table  4:  Summary  of  Characteristics  as  measured  by  the  SllO  and  AFRL  at  300K 


Characteristic 

SllO 

AFRL 

Resistivity  [ohm-cm] 

2.80E-04 

5.60E-05 

Hall  Mobility  rcm^2/(Vs)l 

804 

91.1 

Bulk  Carrier  Cone,  [l/cm'^3] 

2.70E+19 

3.80E+20 

Suggested  Studies 

A  repetition  of  this  experiment  is  merited  sinee  measurements  between  room 
temperature  and  500K  may  still  prove  to  show  surfaee  annealing.  The  following 
suggestions  may  then  be  helpful: 

•  Use  a  probe  station  to  eormeet  samples.  Many  of  the  problems  in  this  experiment 
stem  from  the  deeision  to  pre-attaeh  wires  to  the  samples.  Wires  mean  more 
material  to  be  aetivated,  and  must  be  carefully  chosen.  Additionally,  a  fragile  wire 
and/or  wire  cormection  to  the  contacts  can  break  rendering  an  otherwise  good 
sample  useless.  No  wires  mean  a  smaller  package  to  be  irradiated,  and  less 
cadmium  required  to  shield,  which  is  important  if  the  sample  effects  reactivity  in 
the  reactor. 

•  Allow  ample  time  for  each  measurement.  Measurements  for  this  experiment  took 
about  8  hours  to  complete.  At  a  rate  of  a  sample  a  day,  30  samples  would  take  at 
least  a  month. 

•  Find  a  system  that  is  able  to  complete  both  low  and  high  temperature 
measurements.  Many  systems  only  conduct  high  or  low  temperature 
measurements. 
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Glossary  of  Abbreviations 


AFIT 

Air  Force  Institute  of  Teehnology 

AFRL 

Air  Force  Research  Lab 

CIF 

Central  Irradiation  Facility 

DI 

Deionized 

MOS 

Metal  Oxide  Semiconduetor 

NTD 

Neutron  Transmutation  Doping 

OSUNRR 

Ohio  State  University  Nuclear  Research  Reactor 

PKA 

Primary  Knock-on  Atom 

SiC 

Silicon  carbide 

USAF 

United  States  Air  Force 
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Appendix  A:  Hall  Effect  Procedure 


NIST  posted  the  following  reeommended  procedures  based  on  ASTM  F76  as  a  starting 
point  for  beginning  Hall  Effect  measurements.  [41] 

Hall  Effect  measurements  using  the  Van  der  Pauw  method  requires  an  arbitrarily  shaped 
sample  with  four  contact  on  the  periphery  as  far  apart  as  possible.  A  cloverleaf  geometry 
is  preferred  to  minimize  error,  but  sample  can  also  have  a  square  shape. 

The  following  procedures  correct  for  possible  electrical  non-symmetry  from  non- 
symmetric  contact  placement  by  taking  two  sets  of  measurements,  reversing  the  magnetic 
field  between  sets. 

Part  A)  Van  der  Pauw  measurements  without  Magnetic  Field 

1)  Setup  current  through  sample  with  power  less  than  5  mW  (pref  1  mW) 

2)  Apply  current  Ei  and  measure  V34. 

Note:  I21  refers  to  the  current  entering  contact  2  and  leaving  contact  1. 

V34  refers  to  the  voltage  across  contacts  3  &  4. 

3)  Reverse  the  current  applying  I12  and  measuring  V43. 

4)  Repeat  steps  2  &  3  to  measure  V41,  Vh,  V12,  V21,  V23,  V32 

5)  Calculate  the  Resistance  for  each  I,  V  data  pair  using  ohm’s  law. 

Accuracy  Check:  The  following  relations  should  be  correct  to  within  5% 


D 

_^34 

^21,34 

^21 

R21,34 

=  Ri2,43 

R43,12  — 

R34,21 

R32,  41 

=  R23,14 

Ri4,23  = 

R41,32 

R21,34 

+  Ri2,43  = 

rq 

II 

+  R34,21 

R32,41 

+  R23,14  = 

-  Ri4,23 

+  R41,32 

If  accuracy  check  fails,  check  system  for  possible  sources  of  error. 
Part  B)  Hall  Measurements  with  Magnetic  Field 
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Note:  V24P  refers  to  voltage  measured  in  a  positive  magnetie  field  and  V24N  refers 
to  voltage  measured  in  a  negative  magnetie  field. 

1)  Apply  a  positive  Magnetic  Field 

2)  Apply  current  In  and  measure  V24P 
Apply  current  I31  and  measure  V42P 

Repeat  to  measure  Vnp  &  V31P  using  I42  &  I24  respectively 

3)  Reverse  the  Magnetic  Field 

4)  Repeat  step  2  to  measure  V24N,  V42N,  Vhn,  &  V31N 


Part  C)  Calculations 

Let  the  following  variables  be  defined  for  convenience: 

Vc  =  V24p-V24N  Vd  =  V42P-V42N 

Ve  =  Vi3P-Vi3N  Vf  =  V31P-V31N 


The  sample  carrier  type  can  be  determined  using 
Vc  +  Vd  +  Vg  +  Vp  >  0  :  p-type 
<  0  :  n-type 


Sheet  Carrier  Density 
Ps  = 


Sheet  Resistance, 

Rg  is  found  by  solving  the  equation  =  1 

if  R^  «  Rg  then  Rg  ==  pRjl ,  however  p  is  not  known,  otherwise  it  must 

be  solved  numerically. 

Numeric  Algorithm 

1.  Set  a  relative  error  of  d  =  0.0005 


8-10'MB 

^(Vc  +  Vp)  +  Vp  +  Vp  ) 
8-10'^  IB 

<?(Vc+Vb+Ve+Vp) 
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2.  Calculate  z„  = - 

Jt{R,+R,) 

3.  Calculate  - 

exp(jE._ii?^)  exp(jiz._ii?^) 

4.  Calculate  z,.  =  z._j - - - - 

R^  /  exp(jzz,_ii?^  )  +  RbI  exp(jiz._ii?5 ) 

5.  Repeat  Steps  3  &  4  until  (z.  -  z._^  ^jz.  <  6 

6.  Sheet  Resistance,  R^  =V^i 

Hall  Mobility,  pijj 

or  =  y^Ps^s 

The  following  values  require  the  conducting  layer  thickness,  J  ,  to  be  known 

Resistivity,  p  ■  cm] 
p  =  Rsd 

Bulk  Carrier  Density,  n  or  p  |^carncr/cm^j 
n  =  n^ld  for  n-type  material 
p  =  Psid  for  p-type  material 
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Appendix  B:  Reactor  Time  Calculations 


Appendix  B  steps  through  the  calculation  necessary  to  determine  time  required  to 
receive  a  IMeV  equivalent  fluence  of  lO'"^,  lO'^,  and  lO'^  in  silicon  carbide  in  the  CIF  at 
OSUNRR. 

Heavy  Particle  Collision: 

Heavy  particles  interact  in  much  the  same  manner  as  billiard  balls.  Ignoring  effects  due  to 
Coulombic  and  strong  nuclear  forces,  one  can  easily  derive  an  expression  for  the 
maximum  amount  of  energy  imparted  by  two  particles  of  differing  mass  using 
conservation  of  energy  and  momentum.  The  result  is  the  following  equation,  where  m 
and  M  are  each  particle’s  respective  mass,  and  E  is  the  energy  of  the  incident  particle. 

W.M 

Gmax  =  7 - —2  E  =  (8) 

(m  +  M) 

The  expression  of  the  ratio  of  masses  can  be  expressed  as  simply  a  constant,  where 
convenient. 


The  maximum  energy  imparted  assumes  a  head  on  collision.  As  the  incident  particle 
collides  at  different  angles,  it  will  exchange  less  percentage  of  its  energy.  The  average 
energy  imparted  in  a  collision  between  two  masses  is  considered  to  be  one-half  of  the 
maximum  energy  transferred  in  the  collision,  and  is  modeled  by  equation  9. 


Qavc 


(9) 
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For  a  neutron  colliding  with  silicon  or  carbon,  equation  1  results  in  an  average 

“  0-104 .  Energy  of  30  to  40  eV  is  required  to  displace  an  atom  from  its  lattice 
position  in  SiC.  Neutrons  below  1  keV  rarely  cause  displacements,  thus  rarely  contribute 
to  displacement  dose. 

Monoenergetic  Dose: 

Dose  is  defined  as  the  amount  of  energy  imparted  to  a  system  per  unit  mass.  For  a  fluence 
of  monoenergetic  particles,  dose  is  equivalent  to  that  fluence  {xp)  times  the  scattering 
macroscopic-cross  section  (2^,  probability  of  collision  and  scatter  per  unit  length),  times 
the  amount  of  energy  imparted  per  collision  ( Q).  Since  the  amount  of  energy  per 
collision  changes  with  the  incident  angle,  the  average  energy  per  collision  is  used. 
Macroscopic  cross  section  is  defined  to  be  the  density  of  atoms  times  {N)  times  the 
microscopic  cross  section  (ct^),  resulting  in  equation  3. 

D  =  (10) 

Mulit-Energy  Dose: 

Since  neutron  sources  are  rarely  monoenergetic,  total  dose  is  an  integration  over  energy. 
Since  the  amount  of  energy  imparted  to  the  system  is  Q,  Dose  is  an  integration  over  dQ, 
which  can  be  broken  down  into  a  constant  times  dE .  Thus,  dose  is  reduced  as  follows. 

/1 00 

^^Na,-dE  (11) 

The  OSU  Research  Reactor’s  GIF  facility  has  a  neutron  energy  distribution  such  that 

neutron  flux  between  energies  of  1  eV  and  10  MeV  is  roughly  constant  as  lO'^  - ^ — . 

cm  s 

Since  neutron  energies  below  1  KeV  rarely  contribute  to  the  displacement  dose,  and 
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neutron  population  drops  significantly  at  10  MeV,  the  assumption  is  made  that  only 
neutrons  from  10^  to  10^  eV  contribute  to  dose.  The  assumption  is  also  made  that  N  stays 
constant  since  a  significant  quantity  will  not  transmute  in  the  reactor.  Applying  these 
assumptions,  and  using  flux  (jp)  instead  of  fiuence,  the  dose  rate  (D)  can  be  modeled  by 

D-C_Ni,f°’'Pa,(E)-dE.  (12) 


Thus,  only  the  energy  dependence  of  the  microscopic  scattering  cross  section  of  the 
material  needs  to  be  additionally  known. 


IMeV Equivalent  Fiuence: 


The  1-MeV  equivalent  fiuence  is  the  fiuence  of  1  MeV  neutrons  that  would  contribute  a 
similar  amount  of  dose  as  the  full  spectrum  fiuence.  Assuming  steady-state  conditions,  so 
that  all  variables  are  time  independent,  then  the  1  MeV  equivalent  flux  is  expressed  by 
equation  13. 


^IMeV 


_ ^  ^low _ 

^s@\Mev  '  IT/eV 


(13) 


The  required  time  is  then  represented  by  the  quotient  of  the  target  fiuence  to  the  flux, 
presented  by  equation  14. 

t  =  (14) 

'4’lMeV 


Using  the  microscopic  crossections  from  ENDF  files  for  silicon  and  carbon,  produces 
table  5  as  the  required  times  for  the  desired  fluences. 
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Table  5:  Reactor  time  requirements  for  target  closes 


Fluence 

Power  [kW] 

Time  [s] 

l.OOE+14 

500 

66.67 

l.OOE+15 

500 

666.7 

l.OOE+16 

500 

6667 

Table  5  reflects  the  amount  of  time  required  to  achieve  the  target  fluence  given  the 
reactor  is  run  at  full  power  (500kW).  Since  neutron  flux  scales  linearly  with  power, 
decreasing  the  power  hy  a  given  fraction  increases  the  required  time  hy  the  reciprocal  of 
that  fraction. 
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appeared  minimally  affected  by  neutron  irradiation.  This  suggests  the  creation  of  active  acceptor  defects  decreasing  carrier  concentrations.  N-type  samples 
measured  were  5mm  x  5mm  square  with  Nickel  contacts,  and  irradiated  to  lO'^  and  lO”’  n-cm/cm’of  IMeV  equivalent  neutron  fluence.  Suggestions  for 
continuing  research  include  using  a  probe  station  instead  of  wire  connections  to  samples,  use  a  large  source  current  to  minimize  variance,  and  minimize 
cadmium  shielding  to  reduce  negative  reactivity. 
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